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Abstract We report the discovery of Balmer broad absorption lines (BALs) 
in the quasar LBQS 1206+1052 and present a detailed analysis of the peculiar 
absorption line spectrum. Besides Mg II AA2796, 2803 doublet, BALs are also 
detected in He I* multiplet at AA2946, 3189, 3889 A arising from metastable 
helium 2 3 S level, and in Ha and H/3 from excited hydrogen H I* n = 2 level, 
which are rarely seen in quasar spectra. We identify two components in the BAL 
troughs of Av -~2000 km s _1 width: One component shows an identical profile 
in H I*, He I* and Mg II with its centroid blueshifted by — v c « 726 km s _1 . 
The other component is detected in He I* and Mg II with —v c « 1412 km s _1 . 
We estimate the column densities of H I*, He I*, and Mg II, and compare them 
with possible level population mechanisms. Our results favor the scenario that 
the Balmer BALs originate in a partially ionized region with a column density of 
N-r ~ io 21 ~ 22 cm -2 for an electron density of n e ~ 10 6 ~ 8 cm -3 via Lya resonant 
scattering pumping. The harsh conditions needed may help to explain the rarity of 
Balmer absorption line systems in quasar spectra. With an i-band PSF magnitude 
of 16.50, LBQS 1206+1052 is the brightest Balmer-BAL quasar ever reported. Its 
high brightness and unique spectral properties make LBQS 1206+1052 a promising 
candidate for follow-up high-resolution spectroscopy, multi-band observations, and 
long-term monitoring. 

Key words: galaxies: active — quasar: absorption lines — quasar: individual 
(LBQS 1206+1052) 



1 INTRODUCTION 



Outflows from active galactic nuclei (AGN) play an important role i n galaxy evolution. Ga s 
in AGN may be accelerated to high velocities by therm al pressure (iBegelman et all Il983l ). 
radiation pressur e due to line a nd continuum absorption (jProga &: Kallmanl . 120041 ) . magneto- 
centrifugal forces ( Everettl . 12005). or their combinations. The outflowing gas is potentially crucial 
to maintaining the accretion process by carrying away angular momentum; to chemical evolu- 
tion in AGN by tunneling metal-enriched circum nuclear gas into their hosts ( Hamannl Il999t 
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IWang et~ .,2009): and to regulating the co-growth of the black holes and the bulges of their host 
galax ies by feedback of energy and momentum (e.g. iTremaine et all . 120021 : iFerrarese fe Merrittl . 
l2000h . and t hus quenching the gas supply to the nuclear activity and the host star formation 
process (e.g. iGranato et all . 12004 [dT Matteo et all l2005h . 

AGN outflows often manifest themselves as blueshift ed absorption lines, com monly clas- 
sified into three categories depending on their width Av ( Hamann fe Sabral l2004h : (1) broad 



absorption-lines (BALs) with Av > 2000 km s _1 ; (2) narrow absorption lines (NALs) with 
Av < 500 km s , (3) mini-BALs with Av in between. NALs are ver y common with a de- 
tection rate of about 6 0% in Scyfert galaxies ( Crenshaw et all 1 1999h and 50% in quasars 



(jHamann fe Sabral . 12004ft . BALs are seen in ~ 15% of optically selected quasars^. NALs can orig- 



inate from gas either intrinsic or intervening to AGN, while BALs are generally believed to be 
intrinsic to AGN. All known BAL quasars show high- ionization broad absorption lines (HiBAL), 
including C IV, Si IV and N V. About 15% of BAL quasars present low -ionization broad 
absorption-lines (LoBALs) besides HiBALs, su ch as Al III, Al II, and Mg II (IWevmann et all 
ll99lHReichard et all [2003HZhang et~aTll2010l) . A minority (another 15%) (|Hall et all l2002h of 



LoBAL quasars exhibit additional absorption in Fe II and Fe III from both ground and excited 
levels (FeLoBAL). Even rarer are excited He I* absorption lines, and Balmer absorption lines 
are the extreme rarities. 



Only 5 Balmer-absorption AG Ns have been repo rted thus far: NGC 4151 (jHutch ings et al. 
12002ft . SDSS J083942. 11+38 0526.3 (lAoki et all I20061). SDSS J125942.80+121312.6 (iHalll . 12007" 



SDSSJ102839. 11+450009.4 (IWang et all 12008ft . and SDSS J172341. 10+555340.5 (lAokiboiOl) .' 



Three of the five are formally classified as FeLoBALt@. The high occurrence rate of both Fe 
II and Balmer absorption lines indicates that they are closely related, in spite of the small 
number of such AGN. Their UV spectra are often prof use in excited Fe II* absorption lines 
that are help to constrain the density of absorption gas (|de Kool et al .. 2001). We would have 
a more complete picture of the physical condition and geometry, i ncorporating ioni z ation state 
with He I* as a useful trac er. Interestingly, b oth of NGC 4151 ( Hutchings et all l2002ft and 
SDSS J102839. 11+450009.4 (jWang et all 120081) show He I* NALs, which are the only Balmer- 



absorption AGN with He I* covered by previous spectroscopy with negligible contamination. 
Note that only a handful of He I* absorption AGN have been reported in the literature^. The 
co-occurrence of He I* and Balmer absorption lines is again remarkably high. 

Balmer absorption lines arise from the absorption by n = 2 shell excited hydrogen atoms. 
H I* n = 2 level could be populated by recombination, collisional excitation and Lya resonant 
scattering ((Hail . 120071: fWang et all 12008ft . He I* lines arise from absorption by highly metastable 



level He I* 2 3 S with transition probability A(2 3 S, l 1 ^) = 1.26 x 10~ 4 s" 1 . Transitio ns from 
2 3 S to different upper levels include He I* A10830 in the infrared (jLeighlv et all 1201 ll ). A3889 



in the optical, and AA3189, 2946, 2830, 2764 in the UV. The population of this level comes 
from the recombination of He + ion. A plenty of > 24.6 cV photons are needed to ionize neutral 
helium, while there should not exist too many > 4.8 eV photons to ionize He I* 2 3 S. Diffuse 
stellar background provides too many photons to ionize the He I* 2 3 S' level, and thus these lines 
are not easily seen in the interstellar medium. We c an therefore safely associate the incidence 
of neutral helium line with an AGN environment ( Rudy et all Il985l ). Ignoring second-order 



1 The frequency de pends on th e spec ifi c quasar samp l es and different BAL pick i ng techiques, r angin g 
Tolea et all 12002ft ; iReichard et al.l J2003ft ; iTrump et all J2006ft ; iGibson et all £2009); 



from 10% to 40% ( , 
Gan guly fc Brothertonl £2008)) 



2 The two non-FeLo BALs are NGC 4151 and SDSS J102839.11+450009.4 The forme r show s Fe II NALs 
JHutchings et al.1 . 12002T) . and Fe II NALs are marginally detected in the later JWang et al.1 , 12008ft . 

3 Previo usly known He I* abso r ption AGN inclu de Mrk 231 jBoksenberg et all 119771; iRudv et al.1. Il985f) , 
NGC 4151 JHutchings et al.1 |2002t [Anderson! , 11974ft . Q 2359-124l lArav et all I2001T ). IWang et al.l j200gT) re- 
portcd SDSS J102839. 11+450009. 4 as a remarkable s trong narrow Fe II emission-line quasars with He I* 
AA3889, 3189, 2796 narrow absorption lines. Recently, iLeighlv et al.l J2011ft reported first observation of the 
infrared line He I* A10830, which is 23 times stronger than optical He I* A3889 and could determine a much 
lower column density. 
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effects, reco mbinations of He + to He I* 2 3 S' is mainiy baianced with collisionai de-excitation 
of this level ( Arav et all 1200 lh . Upon that, He I* lines can be used to set a lower limit on the 



He + column density, or a lower limit to H II column density with an assumed abundance. 

Potentially we could have a good chance to explore the properties of AGN Balmer absorption 
line region (ALR). Due to their rarity, however, it is not yet clear whether AGN with Balmer 
ALR form an unusual class, or they are otherwise normal AGN with our line-of-sight (LOS) 
happened to penetrate Balmer ALR, which exists in most if not all AGN, but with an extremely 
small global covering factor (GCF). It is not even known how rare Balmer-absorption AGN are. 
Discovery of more such AGN are obviously helpful for understanding their intriguing nature. 

In this paper, we report the new identification of Balmer BALs in LBQ S 1206+1052, which 
was identified du ring the large bright quasar survey ( Hewett et all [T995h and recently classi- 



fied as a LoBAL dGibson et all l2009h . It is in fact one of ~ 20 Balmer-absorption quasars we 



recently identified duri ng our systematic search in the quasar catalog from the 7th data release 
(ISchneider et all l2010l) of Sloan Digital Sky Survey (SDSS, (jYork et all l2000h L With i band 



PSF magnitude of 16.50, the high brightness of LBQS 1206+1052 makes it an excellent candi- 
date for high resolution spectroscopic follow-up and for long-term variability campaign which 
are crucial for absorption line study. The paper is organized as follows. In Section 2, we analyze 
the SDSS spectrum and measure column densities of H I*, He I*, and Mg II. The results are 
discussed in Section 3. Section 4 is a brief summary of our main conclusions. 

2 SDSS SPECTROSCOPY 

2.1 Continuum and Emission Line Spectrum 

LBQS 1206+1052 was targeted as a quasar candidate for spectroscopy, and was observed by 
SDSS 2.5m telescope on 2003-03-24 with an exposure time of 3165 seconds. It is very bright 
with optical PSF magnitudes of 17.19, 16.63, 16.50, 16.34, 15.92 at it, g, r, i, z bands. The SDSS 
spectrum, retrieved from SD SS Data Release 7 da tabase, is corrected for a Galactic reddening 
of E{B — V) — 0.023 mag (jSchlegel et all Il998| ). and brought to the source rest-frame for 
further analysis with the systemic rcdshift (z = 0.3947 ± 0.0008), re-measured from [O III] 
A5007 as described below. The spectrum is shown in Figure [T] In the rest wavelength range 
of A - 2500 - 6600 A, BALs are evident in Ha, H/?, He I* AA3889, 3189, 2946, and Mg II 
AA2796, 2803 with a velocity range of Aw - - 2000km s -1 . 

Properly modeling of continuum and emission lines is essential to derive absorption line 
parameters. In the observed quasar spectrum, prominent broad and narrow emission lines are 
siting on top of pseudo-continuum, consisting of thermal emission from the accretion disk, 
Balmer continuum from broad line region or disk, the blended Fell lines and high order Balmer 
lines. We fit the pseudo-continuum and emission lines in two steps as follows. In the first step, 
the pseudo-continuum is modeled with three components: a single power law to mimic emission 
from accretion disk, a Balmer continuum plus high order Balmer lines, and blended Fe II 
emission lines. 

Following iDietrich et al.l (|2003l ) and iTsuzuki et al.l ((2006) , the Balmer continuum is ex- 
pressed by F\ = ^3646 x B\(T e )(l — e~ TA ) for A < 3646 A, where i<3646 is the normalization 
coefficient at Balmer edge 3646 A and t\ = tbe(A/Abe) 3 in which tbe is optical depth at 
Balmer edge Abe (3646 A), and B\(T e ) is the Planck f unction at an electron temperature T e , 
which is assumed to be 15, 000 K ([Dietrich et al .. 2003). High order Balmer lines up to n = 50 



are also used in the fit. Relative strengths of these lines ar e fixed using the line emis sivities for 
Case B, T e = 15, 000 K and n e = 10 s cm -3 as calculated bv lStorev fe Hummed (|l995h . and each 



line is assumed a Gaussian profile with a FWHM of 3000 km s -1 . The relative flux of the high 
order Balm e r line s to the flux of the Balmer continuum at the edge is fixed using the results in 
I Wills et al.l ( 19851 ). Thus these blended high order Balmer lines and the Balmer continuum join 



smoothly, forming a pseudo-continuum ( Wills et all 1 1985h . Note that none of the assumptions 
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concerning the Balmcr continuum and high order Balmer lines will have significant effect on 
the measurement of absorption lines because all these features are quite broad. 

The optical and UV Fe II emission are modeled separately with broadened empirical tem- 
plates derived from observed quasar spectra. For optical band long-ward of 3500 A, we use a 
t emplate extracted f r om th e optical spectrum of proto-type narrow line Seyfert 1 galaxy I Zw 
1 IVeron-Cettv et al.l (|2004 ). The template includes both the broad and narrow component, the 
latter also includes transitions from Ni II, Cr II and Ti II. For UV portion short-ward of 3500 A, 
we employ a template built from SDSS J1632+3405 (Zhou et al. 2012 in preparation). SDSS 
J1632+3405 has ultra strong UV Fe II emission, which has a FWHM of 400 km s -1 , only about 
half of that in I Zw 1. The narrow width makes it easy to separate Fe II multiplets from the 
other lines and also among different multiplets. Thus it is a more promising target for building 
UV Fe II template. The relative strength of different UV multiplets ratios are allowed to vary 
during the fit. This is a common technique whe n mod elling quasar Fe II spectra ( Tsuzuki et all 



during tne nt. inis is a common technique whe n mod e 
l2006t iKovacevic et all l2O10t ISameshima et all [20lTh . 

Above continuum components are combined to fit the SDSS spectrum in following windows: 
2855-3010, 362 5-3645, 4170-4260, 4430-4770 , 5080-5550, 6050-6200 A, which are devoid of strong 
emission lines (jVestergaard fc WilkeslkoOlh . The window 3625-3645 is used to constrain Balmcr 
continuum. The best fitted parameters are derived by minimization of x 2 - Line spectrum is then 
acquired by subtracting the best fitted pseudo-continuum from the original spectrum. 

In the second step, emission-lines are measured in this line spectrum. All narrow lines or 
narrow component of broad lines except [O III] and [Ne III] are fitted with a single Gaussian. 
The centroids of the Gaussian are fixed to the vacuum wavelengths of these lines, and their 
widths are tied during the fit. [O III] and [Ne III] show an additional blueshifted component, 
and are fitted with two Gaussians for each of the multiplct. The width and velocity shift between 
the two Gaussians are tied for these lines during the fit, while their normalizations are allowed 
to vary independently. [Ne III] A3869 emission line is affected by He I* A3889 absorption, and 
the pixels affected by He I* A3889 line are masked during the fit according to corresponding 
He I* A3189 absorption profile. Narrow lines included in this fit are [Ne V] A3425, [O II], [Ne 
III] A3869, H(5, H 7 , [O III] A4363, He II A4686, H/3, [O III] AA4959, 5007, [O I] AA6300, 6364, 
Ha, and [N II] AA6548, 6583. 

Broad lines are fitted using as many Gaussians as statistically justified. Broad component 
of H/3 is fitted with 4 Gaussians, while Ha with 3 Gaussians because the red part of Ha line is 
not in the observed spectrum range. Pixels affected by Balmer absorption are masked as done 
for He I* A3889. The fitted profiles of broad Ha and H/3 lines are very similar in the overlapping 
portion in the velocity space, lending credit to our fitting scheme of broad lines. Mg II lines are 
fitted with two Gaussians in the same way as for Balmer lines except that no narrow Gaussian 
is used because there are no explicit narrow Mg II doublet. 

The best fit model is overlaid in Figure [TJ together with individual components. 

2.2 Absorption Line Spectrum and Column Densities 

We use the best fit model to normalize the absorption line spectrum. We first subtract model 
narrow emission lines from the observed spectrum before normalization, since absorption gas 
does not cover the narrow emission line region (NLR) in all well studied BAL quasars. Mg 
II doublet and Balmer absorption lines are superposed on top of broad emission lines. They 
should be normalized in different ways d epending on whethe r the absorbing gas covers the 
broad emission line regions (BLR) or not ( de Kool et all l200ll) : (1) The observed data should 



be divided by the sum of the continuum and broad emission lines if it does. (2) If not, the broad 
emission lines should be subtracted before the observed data normalized by the continuum. We 
adopt the second normalization scheme for the following hints. We bloat Figure 1 and display 
the close-ups in Figure 2 focusing on the two Ha and Mg II regimes. It can be seen clearly 
that the modeled broad emission lines just brush the tip of the BAL trough in both of Ha and 
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Mg II lines. The fact can be most naturally interpreted in the second scheme with an optical 
depth t > 1 at the centroid of both Ha and Mg II BAL troughs. Otherwise we would need 
to have two coincidences at once, which is very unlikely if not impossible. The interpretation 
is also consistent with the marginal detection of Ca II K in absorption, which suggests a large 
optical depth in Mg II BAL doublet for any reasonable abundance. 

The normalized absorption line spectrum is shown in Figure [31 and the apparent optical 
depth profile of Mg II AA2796, 2803 doublets is shown in Figure [U The normalized flux at the 
line centroid approaches zero for Mg II and Ha lines. Note that Balmer emission line profiles 
at the absorption line position is somewhat uncertain due to lack of independent constraint. 
We use 1-er error of normalized flux to deduce a lower limit for apparent optical depth of 
the deepest pixels. The Mg II profile apparently show two components, one broad component 
with corresponding Balmer and helium absorption lines. The other one is relatively narrow and 
shows He I* lines, possible Ca II A3934 absorption line and little Balmer absorption line. We 
simultaneously fit the apparent optical depth t(v) of all absorption lines with two Lorentzians 
for each line, which is similar to Voigt profiles for small b values. The line centroid of other lines 
are tied to that of Mg II lines and the line width is assumed to be equal for the same component 
during the fit. The best fitted model for Mg II is also shown in Figure SI and overlaid in Figure 
|3J The fitted EWs for each line are shown in Table [T] The low velocity component is centred at 
— 726 km s" 1 and with a FWHM of 603 km s , it is referred to as component 1 (CI) hereafter, 
while a high velocity component with a line centroid of —1412 km s _1 and a FWHM of 208 km 
s" 1 are referred to as component 2 (C2) hereafter. 

The CI component shows similar profiles in Mg II, He I* and Balmer absorption lines. 
There is only one pixel in Ha line that is inconsistent with our model, which can be attributed 
to imperfect emission line modeling. The r ratio of 2796 to 2803 for CI is 2:1, which is consistent 
with that the absorption gas fully cover the continuum source. To check this, we compare the 
ratios of Balmer lines with their theoretical values. From Table [I] we can see that the Balmer 
absorption line EW ratio EW(Ha) : EW(H/3) = 8.19^;-^ is consistent with theoretical /A ratio 
7.26, thus we favour that the absorbing gas of CI fully covers the continuum source and the 
lines are not saturated. Under this assumption, we d erive column de nsities of various ions from 
their equivalent widths W\, using the equation from I Jenkinsl ( 19861 ): 



We get a column density for Mg II from Mg II A2803 line: N Mg+ = (1.41 ± 0.07) x 10 14 
cm" 2 , for hydrogen at n = 2 from H/3 line: N H (n = 2) = (1.37 ± 0.10) x 10 14 cm" 2 and for 
helium from He I* A3889 line N Hc i 2*s = (8-99 ± 0.50) x 10 14 cm" 2 

Component C2 shows consistent absorption trough in Mg II and He I*. A weak line near 
3930 A may be attributed to the Ca II K absorption line of this component (blueshifted by 100 
km s _1 relative to other lines). The Mg II absorption line is seriously saturated as indicated 
by the doublet ratio of Mg II, EW2796/EW2803 ~ 1- Thus we can only get a lower limit of 
Mg II column density using equation ([!]), A^ Mg + > (3.43 ± 0.34) x 10 14 cm~ 2 . Alternatively, 

we may integrate the apparent optical depth profile N = ^ifjjz J t v dv and get a column 
density of N Mg + = (5.84 ± 0.45) x 10 15 cm -2 , which is likely close to the real value. The 
Ha line of C2 is not detected significantly. We derive an upper limit on its column density 
Nno(n = 2) < 9 x 10 12 cm" 2 . The He I* lines in this component are generally weak, we get a 
column density of iV Hc 1 2 3 s = (1-25 ± 0.22) 

3 DISCUSSION 



The overall rarity of Balmer absorption line quasar indicates that either the conditions for 
outflowing gas to maintain a sufficient population of n — 2 level hydrogen atoms must be very 
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harsh. We will discuss the conditions required to produce the Balmcr absorption lines in CI 
with different mechanisms. 

If the n = 2 level hydrogen atoms of this component are caused by collisional excitation, 
they must exist in H I region or probably the partially ionized region. The observed Mg II 
absorption lines suggest that there is partially ionized absorption line region. Mg + ions are 
created by photons with energy greater than 7.6 eV and destroyed by photons with energy 
greater than 15.0 eV, thus they mainly survive in outer part of H II region and mostly in partially 
ionized region. Assuming that Mg + is the dominant species of Mg and a solar abundance with 
[Mg/H]=-4.4, we got a lower limit on the total hydrogen column density of 3.54 x 10 18 cm' 2 for 
the partially ionized zone. If thermal equilibrium is achieved, then using Boltzmann equation: 
nilri\ = 4exp(— 10.2eV/fcT), in which n\ and n 2 are densities of n = 1 and n = 2 level hydrogen, 
one yields readily a column density of n = 2 hydrogen of 6.8 x 10 13 cm -2 , which is close to 
the observed hydrogen column density (1.37 ± 0.10) x 10 14 cm" 2 , at a typical temperature 
T~ 10 4 K for photo-ionized gas. However, thermal equilibrium requires that collisions with 
thermal particles are the dominant process for excitation and de-excitation. This yields a lower 
density limit n c = A 2 i/[(72i(l + tqi)] = 8.7 x 10 16 /(1 + t ;) cm" 3 accounting for the resonant 
scattering of Lya photons, where tqi is the optical depth of Lya from the central of the region. 
For a gaussian distribution turbulence, we can write To; ~ 7.6 x 10 6 feg 1 A^h<i,22 ~ 10 3 , where 
63 ~ 1 is the turbulent velocity of gas in units of 10 3 km s _1 and Aho 22 neutral hydrogen 
column density in 10 22 cm" 2 . This gives a density too high for any reasonable line absorbing 
gas. 

When density is lower than the critical density discussed above, detailed equilibrium of 
n = 2 level must be considered. Transition of 2s state to the ground level is forbidden with 
a small Einstein coefficient A = 8.23 s _1 via 2-photon process, while transition from 2p to Is 
is allowed. The two levels have very different excitation processes, therefore, we will consider 
them separately in the following analysis. 

Let us first consider the equilibrium of 2p level. In the case that absorption gas is opti- 
cally thin to Lya, the equilibrium between the recombination to n = 2 and radiative decay 
from 2p to 2s yields n 2p /riH+ = Q. e sn e /-Mp,2^ However, the gas is likely very optically thick 
to Lya. In this case, each L ya photon is scattered 1 + To; times on average before escaping 
( Osterbrock fc Ferlandl . [20061) . This leads to the 2p population larger by a factor of 1 + T01 than 

7.3 x 10"V,6^h",22^3 1 - To 



above optically thin formula, i.e., n2 P /n H + ~ To\a c sn e /A2 P .2s 
get the observed n = 2 column density, we require 



A, 



12 ^ 1.4/ H+ (1 - / H +)" ~'~0 3 ' n H,6 



l-W/a -1/2 



(A n=2 /1.4 x 10 1 



2) l/2 5 



(2) 



where / H + is fraction of ionized hydrogen. As discussed in Section2.2, it is very likely that the 
BALR in LBQS 1206+1052 fully covers the accretion disk but does not cover the BLR. This 
suggest 10 14 cm ~ r A D < ^balr < ^blr ~ 10 17 cm, where r A D, ^balr and r B LR arc the size 
of accretion disk, broad absorption line region, and broad emission line region, respectively. 
We infer a hydrogen column density of Ah ~ 10 21 " 22 cm" 2 of the partially ionized gas for an 
electron density of nn ~ 10 6 " 8 cm" 3 . For a higher ijh and lower Ah, the BALR would not 
be large to fully cover the accretion disk, while the BALR would be large enough to cover the 
BLR for a lower n.H and higher Ah- The inferred Ah value is much larger than the thermal 
equilibrium value. 

Next, we consider population in the 2s level. The equilibrium equation of 2s level can be 
written as 



n e nii+a e B^s + "H« e q , i s -2 S 



n 2s 



Mi + n e (q 2s -2 P + <72.s 



-nil 



a y L u 
4:Ttr 2 hi> 



dp 



(3) 



4 As wc discuss next, the absorbing gas is probably in partially ionized region, where the typical temperature 
gnificantly lower than '. 
excitation of 2p level here. 



is significantly lower than 10 4 K, and qi s -2p oc exp(— Tjr) <S ct c ff in this region, thus we ignore the collisional 
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in which we ignored the effect of 2p level for the time bcin^£]. The first term on the left hand 
represents recombination to 2s level, while the second term represents collisional excitation. 
Given that a e ff,2s and q\ s -2s are in the same order of magnitude (a e e,2s = 0.838 x 10' 13 cm' 3 
s' 1 , <7i s -2s = 1-67 x 10' 13 cm' 3 s' 1 at T=10,000 K), the ratio of the first to second term 
is mainly determined by the ratio of njj+ to n^. In H II region where riH+/ n H is high, the 
recombination mechanism dominates the population of n = 2 level, while in H I region or PIZ 
where n H + /njj is low, the collisional excitation dominates the population of n = 2 level, we 
treat them separately. 

In H I or PIZ region, we ignore the first term on the left hand of equation The third 
term on the right hand is the photoionization of n = 2 level, in which isi is the frequency at 
Balmer edge A = 3646 A. It can be shown that the third term is negligible by substituting the 
integration term with U = \ v L v / '(Airr 2 '■ cnnhv)dv and isq = Av\ and assuming a power low index 
of optical-UV continuum a v = —0.44 (jVanden Berk et "ail . 120011 1. If density n e of absorption gas 
is sufficiently low compared to critical density n c = A2S-+1S /q2S-+2P ~ 1.5 x 10 cm' 3 at T = 
10,000 K (Osterbrockfe Ferland 2006), the equilibrium gives n, 2 s/riH = n e qis^2S I 'A-is-txs < 
<?is^2S /?2S->-2P) otherwise if the density is moderately high, the second route will dominate, 
which gives n2s/ n e = qis^2S /q2S^2P, which is around 3 x 10' 10 at a temperature of 10, 000 K 
0. This requires a total hydrogen column greater than 5 x 10 23 cm -2 to generate the observed 
Balmer absorption lines EW in Q 1206+1052. Clearly, collisional excitation is mush less efficient 
in populating n = 2 level than Lya resonant scattering in partially ionized gas in AGNs if gas 
density is high. 

Alternatively the n = 2 level population could arise from ionized gas via recombination, 
and we ignore the second term on the left hand of equation ((3|) . As in the collisional excitation 
dominated case, we further distinguish cases when densities are below or above the critical 
density 10 4 cm -3 set by collisional angular-momentum transition from 2s to 2p state. If the 
density is higher than that, the 2s state is populated by recombination and depopulated via 
collisional excitation to 2p state, which radiativcly decay to n = 1 level. The equilibrium gives 
"-2s/«h+ ~ a 2s/q2s->2 P ~ 0.44 x 10' 10 at T = 10,000 K. Otherwise if the density is lower 
than that critical density, the equilibrium is built through the balance between the two-photon 
continuum emission of 2s — > Is and the recombination to 2s state, which gives n2 S /nH+ ~ 
n e Oi2s l 'A2s-+is < <^2s I '(Z2s-»2p ~ 0.44 x 10~ 10 . The last equation in above estimations utilizes 
the fact that when this route dominates the equilibrium, n e < n c = A2 S ->i s / '<?2s->2p- Thus if 
indeed the level two population is from the HII zone, the observed n = 2 column require a total 
hydrogen column of at least 3.1 x 10 24 cm' 2 . The column density appears too large for the 
absorbing gas. 

To summarize, the resonant scattering of Lya is much more efficient in populating n = 2 
level in the partially ionized region if gas density is high. In this case, the column density 
required to account for the Balmer absorption lines in CI are in the range of 10 21 to 10 22 cm' 2 
for gas density between 10 6 to 10 8 cm -3 and decreases with square roots of the gas density. At 
much low densities, collisional excitation may take over the resonant scattering pumping, and 
it requires a column density at least 5xl0 23 cm' 2 . In the H II region, recombination can lead 
to the required n = 2 level at very high column density Nn+ = 5 x 10 24 cm' 2 . 

A thick H II region is not consistent with the observed helium absorption lines. He I* 
absorption a rises from absorp tion by metastable He I* 2 3 5* state. The equilibrium equation of 
this level is (|Arav et all 1200 ll) : 



nn c+ n e a T = n 2 3 S 



/CO r 



(4) 



5 Had 2p been significantly populated due to Lya trapping in dense gas as discussed in last paragraph, 
because of large collision strength, 2p and 2s level will reach thermal equilibrium, that is rt2 s : «2p = 1:3. 

6 IWang et alJ {2008) used an incorrect g2S->2P value resulting an incorrect value n2s/riH ~ 0.02 
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where, vq is the threshold frequency at ionization energy 4.77eV of 2 3 S state, and q2 3 s,2 1 s, 
Q2 3 S2 1 Pi 9ci are corresponding collisional rates. The left hand term represent recombination 
of He II ions with electrons. On the right hand of the equation, the first term is radiative 
transition from 2 3 S to l 1 ^ , which is a forbidden transition and can be neglected; the second 
term are collisional de-excitation of 2 3 S I to 2 1 S and 2 1 P; and the third and forth terms represent 
collisional and radiative ionization of 2 3 S respectively. Neglecting terms other than collisional 
de-excitation, we got 

n 2 3 S (He°) = 5.8 x lQ- 6 ^ 119 

n He + i + siioi;- - 51 ^ 1 

where T 4 = T/10 4 K ~ 1 for photoionizcd gas. At a reasonable gas density, the observed He 
I* 2 3 S column density implies a column density of H II zone of only 10 21 cm" 2 . This column 
density is much lower than required for excitation of n = 2 level via recombination, but is 
consistent with Lya pumping as the main excitation mechanism for n = 2 population. 

A thick H I region can produce too larger column density of Mg II. It is shown that the 
fraction of Mg + in the PIZ of AGN ionized region ranges from 10 -3 cm -3 close to H II-H I 
transition zone to near dominated species deep inside the H I region. To meet the observed Mg 
II column density, the PIZ zone can not be thicker than a few 10 21 cm -2 , which rules out the 
possibility of collisional excitation without Lya pumping. Therefore, Lya pumping remains the 
only possible excitation mechanism at least for the Balmer absorption in this object, and the 
density is likely higher than 10 8 cm" 3 . 

For C2, Mg II lines arc seriously saturated. The column densities derived from integrated 
Mg II apparent optical depth profile set a lower limit on hydrogen column density of Ah = 
6.23 x 10 19 cm -2 . On the other hand, a lower limit of the column density (Ah ~ 10 21 cm" 2 ) 
can be derived from the He I* absorption lines with an analysis similar to that for C2. Although 
Mg II column density is much larger than that of CI and He I* column density is similar, the 
Balmer absorption is not detected with an upper limit of 9 x 10 12 cm" 2 . 

As discussed above, n — 2 hydrogen is populated much more efficiently via Lya pumping 
process than with other routes. The process works well in relatively dense material and in 
partially ionized medium. Lack of the Balmer absorption lines can be attributed to any of the 
following cases ([2]): the total column density of partially ionized zone is too low; the ionization 
of the gas ma y be too large; t he density of gas is too low, or any combination of the three. As 
pointed out in lLu et ahl ( 2008h . a trace amount of Mg + may originate from an ionized region with 



log U > —1.5 where Al 2+ and C 3+ live. A total hydrogen column density of Ah ~ 10 22 cm" 2 is 
required to generate the observed Mg + if this is the case, although the column is much larger 
than the lower limit set by helium absorption, we can not rule out this possibility for now. 



4 SUMMARY AND FUTURE PROSPECTIVES 



We present a detailed analysis of absorption line systems in the quasar LBQS 1206+1052, the 
brightest quasar showing Balmer absorption lines Ha and H/3. From its medium resolution SDSS 
spectra, we found that there are mainly two absorption components. Component 1 is blueshiftcd 
by —726 km s" 1 and with a FWHM of 603 km s" 1 , showing corresponding absorption in 
Balmer, He I* and Mg II. This component has full coverage of the continuum source but not 
BLR. Component 2 is blueshifted by -1412 km s" 1 and with a FWHM of 207 km s" 1 . The 
Mg II lines are seriously saturated for this component, while helium lines are generally too 
weak to facilitate partial coverage analysis and Balmer absorption lines are not detected. If 
the H(n = 2) is the outcome of recombination, a typical ionized hydrogen column density of 
several 10 24 cm" 2 is needed to generate the observed absorption in Balmer series. Collisional 
excitation of 2s level without Lya resonance scattering in play also requires an column as high 
as 5 x 10 23 cm -2 . These column density are inconsistent with either Mg II or He I* absorption 
lines. Thus, Lya scattering in a partially ionized region is the main mechanism of populating 
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n = 2 level hydrogen, and the required column is 10 21 to 10 22 cm -2 for gas density between 10 6 
to 10 8 cm -3 and decreases with square roots of the gas density. These harsh conditions required, 
i.e. the large column density and density, may help to explain the overall rare p opulation of such 
objec ts, only ~ 20 (Ji et al. 2011 in preparation) out of 105783 DR7 quasars ( Schneider et all 
120101) . 

Due to the medium S/N and resolution of SDSS spectra, we cannot restrain further about 
the nature of the absorber. However, the bright magnitudes make LBQS 1206+1052 the most 
promising candidate for high resolution spectroscopic observation, which will provide resolved 
velocity profile for pixel-to-pixel based covering factor and velocity structure analysis. The 
Earth's atmosphere cutoff wavelength 3200 A corresponds to a rest-frame wavelength 2285 A at 
z = 0.4 of the quasar. Given the high rate of detecting Balmer absorption in FcLoBALs, we 
expect to detect absorption lines in UV1 multiplet arising from excited levels of Fe II ions 
from ground-based faci lities, such as Keck a nd MMT, which may help to determine the density 
of the absorption gas ( Korista et all 120081) . Long-term monitor is also helpful to explore its 
physical/geometrical properties. 
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Fig. 1 SDSS spectrum of LBQS 1206+1052. Prominent absorption features are 
marked as vertical dashed line and labeled on the top. Model fits to the spectra 
are shown as green line. Different components of the spectral model are also shown. 




Fig. 2 Close-ups of Figure 1 bloated in the Ha (upper) and Mg II (lower) regimes. 
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Fig. 3 Normalized intensity assuming that absorption gas cover only the continuum 
source but not the BLR. From top to bottom are Mg II A 2803, Balmer absorption 
lines and He I* absorption lines, Ca II A3934. Different lines arc plotted using different 
colors as indicated by legends on the lower right of each panel. Best fitted models of 
absorption profile are shown as green lines, only fits for Mg II, Ha, He I* A3889 and 
Ca II A3934 are shown for clarity. The two components are shown in red dot-dashed 
lines for Mg II fit, whose centres are marked as black vertical lines.. The absorption 
model is oversimplified as judged by the central pixels of Ha line, yet has little effect 
on our analysis. 
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Fig. 4 Apparent optical depth for Mg II AA2796, 2803 doublet as function of velocity 
relative to rest-frame 2803 A. The pixels marked as up-arrow are lower limits of optical 
depth values. These pixels has normalized flux approaching zero and 1-er errors of 
normalized flux are used to get the optical depth limit. Various components are also 
plotted using different lines and colors as indicated by legends. 
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Table 1 Absorption line parameters 



Line 


Wavelength(A) 








Ha 


6564.41 


0.6400 


27.88 ±2.16 


< 0.70 


H/3 


4862.68 


0.1190 


3.41 ±0.25 


0.15 ±0.08 


H7 


4341.68 


0.0446 


0.46 ±0.13 


none 


He I* 3889 


3889.75 


0.0278 


3.35 ±0.19 


0.47 ± 0.08 


He I* 3189 


3188.66 


0.0110 


1.46 ±0.13 


0.19 ±0.06 


He I* 2946 


2945.96 


0.0053 


0.91 ±0.13 


0.15 ±0.06 


Mg II 2796 


2796.35 


0.6160 


6.00 ±0.31 


7.29 ± 0.73 


Mg II 2803 


2803.53 


0.3060 


3.01 ±0.15 


7.31 ±0.73 


Ca II 3934 


3934.78 


0.6270 


none 


< 0.24 



" Oscillator strength are from NIST Atomic Spectra Database (http:/ /physics. nist.gov/PliysRcfData/ASD/ ). 
b Equivalent widths for component 1, which is a Gaussian with centroid at -726 km s -1 and FWHM 603 km 



c Equivalent widths for component 2, which is a Gaussian with centroid at -1412 km s 1 and FWHM 208 
km s^ 1 
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